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Abstract. The field-induced reorientation of the magnetization of ferromagnetic films is treated within the 
framework of many-body Green's function theory by considering all components of the magnetization. We 
present a new method for the calculation of expectation values in terms of the eigenvalues and eigenvectors 
of the equations of motion matrix for the set of Green's functions. This formulation allows a straightforward 
extension of the monolayer case to thin films with many layers and for arbitrary spin and moreover provides 
a practicable procedure for numerical computation. The model Hamiltonian includes a Heisenberg term, 
an external magnetic field, a second-order uniaxial single-ion anisotropy, and the magnetic dipole-dipole 
coupling. We utilize the Tyablikov (RPA) decoupling for the exchange interaction terms and the Anderson- 
Callen decoupling for the anisotropy terms. The dipole coupling is treated in the mean-field approximation, 
a procedure which we demonstrate to be a sufficiently good approximation for realistic coupling strengths. 
We apply the new method to monolayers with spin S > 1 and to multilayer systems with S = 1. We compare 
some of our results to those where mean-field theory (MFT) is applied to all interactions, pointing out 
some significant differences. 
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1 Introduction 

In this paper, we extend our earlier investigations KM 
on the reorientation of the magnetization of a ferromag- 
netic monolayer with spin 5 = 1 to multilayer systems 
with arbitrary spin, S. The components of the magneti- 
zation as functions of temperature and film thickness are 
calculated within the framework of a many-body Green's 
function theory, allowing the direct calculation of the mag- 
netic orientation. Furthermore, we derive and apply a non- 
perturbative expression for the temperature dependence of 
the second-order single-ion anisotropy by minimizing the 
free energy with respect to this orientation angle. 

For convenience, mean-field theory (MFT) is often ap- 
plied to such problems, either by diagonalisation of a single- 
particle Hamiltonian H, or by a thermodynamic pertur- 
bation theory 0. We emphasize that this approximation 
completely neglects collective excitations (spin waves) 
which are known to have a much greater influence on 
the magnetic properties of 2D systems than on 3D bulk 
properties. In fact, recent calculations on a trilayer sys- 
tem p demonstrate that MFT is incapable of accounting 
for the induced magnetization observed in coupled lay- 
ers unless an unrealistically large interlayer coupling is 
postulated. On the other hand, the many-body Green's 



function techniques, which take the collective excitations 
approximately into account, can explain the experimen- 
tal observations assuming an interlayer exchange coupling 
of reasonable size. We point out that the dependence of 
the magnetization and the Curie temperature on the film 
thickness has also been studied within Green's function 
theory g. Also, in references and |^ Green's function 
techniques are applied to magnetization problems. How- 
ever, in all these references only a single component of 
the magnetization is considered. Therefore, a reorienta- 
tion of the magnetization as a function of temperature, 
film thickness, or magnetic field cannot be calculated. A 
reorientation of the magnetization is considered in refer- 
ence p, where Green's functions are used after a Holstein- 
Primakoff mapping of the spin operators to bosons. In or- 
der to solve the non-Hermitian eigenvalue problem, right 
and left eigenvectors are applied for its solution, sinrilar to 
what is done in the present paper. The theory is, however, 
only valid at low temperatures. Another method for the 
treatnrent of the magnetic reorientation for all tempera- 
tures is a Schwingcr-Boson theory |1^, as an alternative 
to the Green's function method of this paper. 

In the present work, we treat the field-induced reori- 
entation of the magnetization for all temperatures of in- 
terest. Since expectation values of all three components 
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of the spin operator are considered, a corresponding set 
of Green's functions must be defined. We introduce a new 
method for the calculation of the expectation values, which 
utilizes not only the eigenvalues but also the eigenvectors 
of the (non-symmetric) matrix governing the equations of 
motion for the Green's functions. This formulation is more 
compact than the usual one B and, most importantly, 
suggests a practicable way of treating the multilayer case 
for arbitrary spin. We make no attempt to go beyond the 
Tyablikov (Random Phase Approximation: RPA) decou- 
pling for the exchange terms, since we have shown in ref- 
erence PI that this decoupling scheme for a monolayer 
with spin S = 1/2 compares well indeed with an 'exact' 
Quantum Monte Carlo (QMC) calculation. The single-ion 
anisotropy term is decoupled with the Anderson-Callen 
method |13|, because, as shown in 0, other single- ion 
decoupling schemes, e.g. that of Lines Q, lead to diffi- 
culties when calculating the magnetic reorientation. Fur- 
thermore, we include the magnetic dipolc coupling, which 
is treated within a simplified (non-dispersive) approxima- 
tion, which corresponds to its mean field treatment. The 
dipole coupling was not considered in our earlier work Q , 
which was restricted to the case of a monolayer with spin 
5'=1. 

The paper is organized as follows. In Section 2 the 
Green's function formalism is outlined. The eigenvector 
method is then presented for the monolayer case and it 
is shown that this leads to a transparent extension to the 
case of many layers. Section 3 deals with the results. First 
the effect of the dipole coupling on the magnetic reorien- 
tation is discussed using the monolayer with spin 5 = 1 as 
an example. Secondly, we discuss monolayers with spins 
S > 1. Thirdly the formalism is applied to the case of 
many layers. Finally Section 4 contains a discussion of 
the results and an outlook for further investigations. In 
Appendix A, different approximations for the magnetic 
dipole coupling are investigated. Details of the formalism 
for 5 > 1 are derived in Appendix B. 



2 The Green's function formalism 

In order to study the field-induced magnetic reorientation 
of a ferromagnetic thin film we investigate a spin Hamil- 
tonian consisting of an isotropic Heisenberg exchange in- 
teraction, Jki, between nearest neighbour lattice sites, a 
second-order single-ion lattice anisotropy with strength, 
K2,k, the magnetic dipole coupling with strength, g^^, and 
an external magnetic field, B = (S^, B^, B^), 

^=-^E^«S,.S,-5]if2,fe(5f)2 

<kl> k 



Y.(\B-St + \B^S^+B^-Si 



1 \r^ 9ki 



.^Y.-^[rliSk-Si-'i{Sk-Vki){Si-Vki) 
^ ki '^ki 

Here the notation Sf = Sf ± iSf and B^ = B"" ± iBV 
is introduced, k and I are lattice site indices, and {kl) indi- 



cates summation over nearest neighbours only. Each layer 
is assumed to be ferromagnetically ordered (coUinear mag- 
netization) , whereas the magnetization of different layers 
need not to be coUinearly aligned. Furthermore, inhomoge- 
neous systems can be considered which are characterized 
by different layer-dependent coupling constants and mag- 
netic moments. We do not include a fourth order uniaxial 
anisotropy term, —'}2,k-^kA{Sf.Y, because it is difficult 
to find a proper decoupling of this term in the equations 
of motion for the Green's functions. This means that the 
formalism of the present paper is only adequate for the 
physical situation in which such a term is of no impor- 
tance. 

As in reference g, we introduce the set of thermal 
Green's functions in the spectral representation 



G-7(c.) = ((5f;(5|)™(5-)")). 



(2) 



where lo denotes the energy, and t] = ±1 refers to the com- 
mutator (j] = —1) or anti-commutator [rj — +1) Green's 
functions, respectively; n > \ and ?7i > are integers, i 
and i denote lattice sites. In order to obtain a closed set of 
equations of motion for the Green's functions, we treat the 
exchange term by a generalized Tyablikov (RPA) ||lj| de- 
coupling, and the anisotropy term by the Anderson-Callen 



decoupling [12 



In our previous work the corresponding thermal corre- 
lation functions 



c""'" = {{s'nsjTSf) 



(3) 



have been obtained by applying the spectral theorem [n4| . 
Because of vanishing eigenvalues it is important to use 
the spectral theorem including the term obtained from 
the anti-commutator Green's functions 






(4) 



Together with the so-called regularity conditions, which 
are derived from the fact that the spectral representation 
of the commutator Green's function must be regular for 
a; = 0, we have derived a set of coupled equations for 
the correlation functions. The solution yields the compo- 
nents of the magnetization, thus determining directly the 
reorientation angle of the magnetization induced by the 
applied external field. In the present paper, we rederive 
these equations by a method which utilizes the eigenvec- 
tors as well as the eigenvalues of the matrix determining 
the Green's functions. This more compact formulation fur- 
nishes a practicable way to treat the multilayer case and 
general spin quantum numbers S. It is didactically ad- 
vantageous to demonstrate this new method first for a 
monolayer; the generalization to the multilayer case then 
follows in a straightforward and transparent way. 



2.1 The eigenvector method for the monolayer 

The equations of motion for the Green's functions in the 
spectral representation read 

u g:;.'(7(^) = <y(7 + {{[STM-i: {s^.nsjT)) , (5) 
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with the inhomogeneities 

= {sT{sir{sjr + visinsjrsf) , (6) 

where (••■)= Tr(- • • e"^'^) with (3 = l/fc^T and ks Boltz- 
mann's constant, and rj = +1 or —1, respectively. 

The higher Green's functions in equation (|^) due to the 
exchange interaction term are decoupled by a generalized 
Tyablikov (RPA) decoupling Ol 



rf a Cf/3 , / QZ \ni f Q—\n 



/3,mn |^/q/3\ y-ia.mri 



{{SrS^,;iS]riSrr)) ^ (S?) G^^f " + (^0 Gj"" . (7) 

In reference H] the proper inclusion of the single-ion aniso- 
tropy with Green's function techniques was thoroughly 
discussed in connection with the magnetic reorientation. 
Accordingly, we choose the Anderson-Callen |l3,p| decou- 
pling for the treatment of the anisotropy terms: 



{{S^ S'j + S^ Sj_ )) 



2{Sn{l-^[S{S + l)-{StS^)])G. 



±,T?in 



(8) 



Because we are interested in laterally periodic sys- 
tems we perform a Fourier transformation to the two- 
dimensional wave vector space k. Introducing vectors for 
the Green's functions, G™"(k, a;), and for the inhomo- 
geneities, A™", 



ran 



/G+'™"(k,c.) 
G7'(k,c.)= G--'""(k,w) I , A™"= I A-. 

\G^'™"(k,w)y \A^^™" 

the equations of motion, which are derived in detail in 
reference fei , can be written in a compact form 



(cji-r)G™" = A™, 



(10) 



il^ 





-H+ 





~H'' 


H- 


\H- 


2-" 






where 1 is the unit matrix and the non- symmetric matrix 
r is given by 



(11) 



with the abbreviations 

H"=i?" + (5")J(<Z-7k), a = +,-,z, 

Z = B^ + A'2<Z>, (12) 

For a square lattice with a lattice constant taken to be 
unity, one obtains 7k = 2 (cos kx + cos ky), and g = 4 is the 
number of nearest neighbours. Note that A™" = A™"(k) 
depends on the wave vector k for rj = 1 but not for 77 = —1. 
Similar to the exchange coupling, the long-range dipole 
coupling also induces a momentum dependence into the 



magnon dispersion relation (jj(k). Due to the oscillating 
lattice sums the consideration of the k- dependence of 
this coupling obtained e.g. in RPA is fairly complicated 
and time consuming. Therefore we accept for the present 
calculations an approximate description of the dipole cou- 
pling in the dispersion relation, in particular its k- de- 
pendent terms are neglected, which is equivalent to its 
mean field approximation. In Appendix A, we show that 
for dipole interactions small compared to the exchange 
coupling, which is the case for the ferromagnetic 3d- tran- 
sition metals, this approximation is satisfactory. The ap- 
proximation of Appendix A merely leads to a renormal- 
ization of the external magnetic field components B^ and 
B^ , which for the ith atomic layer in the case of a multi- 
layer with N layers reads 

N 



Bf ^ B^ +Y^ g,, {Sf) T\^'^\ 

N 



(13) 



i=i 



where the lattice sums for a two-dimensional square lattice 
are given by (n = \i ~ j\) 



^" = E 



Im 



{P + ni^ + n2)5/2 



(14) 



The indices {Im) run over all sites of the square j-th layer, 



excluding the terms with I 



{^ — Q. For the mono- 



layer {N — 1) one has i — j, and one obtains in partic- 
ular T" ~ 4.5165. As can be seen from equations (Esl), 
the dipole coupling reduces the effect of the external field 
component in ^-direction and enhances the effect of the 
transversal field components B . 

We now introduce a transformation which diagonalizes 
the matrix T 



LrR = o 












(15) 



where the eigenvalues turn out to be loq = 0, lu± == iii'k 
with E^ — H^H^ + H^H^-. The transformation matrix 
R and its inverse R^^ = L are obtained from the right 
eigenvectors of F as columns and from the left eigenvectors 
as rows, respectively. These matrices are normalised to 
unity: LR = RL = 1. Note that due to the non- symmetric 
matrix F, equation (ll]), one has in general R^^ ^ R'^, 
R"^ being the transposed matrix. 

For the monolayer, the transformation matrices can 
be constructed analytically; the right eigenvectors are ar- 
ranged so that the columns 1,2, and 3 correspond to the 
eigenvalues 0, +-Ek, and — i?k, respectively: 



R 



H+/H' 


-{H^+E^)IH- 


-{H--E^)/H 


H- /H'' 


-{m ~ E^)/H- 


-{H-+E^)/H 


1 


1 


1 



(16) 
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/ 

-4 


2H-H^ 


2H+H^ 


4:H^H'' 




-H+{E^ + H^) 


2H-H+ 
2H-H+ 



(17) 



Similarly, the left eigenvectors are arranged so that rows 1, 
2, and 3 correspond to the eigenvalues 0, +£'k, and — £'k, 
see equation (^ above. 

Multiplying the equation of motion (HG) from the left 
by L and inserting 1 = RL one finds 

(tjl-n)LG;"" = LA™" . (18) 

Defining ^^"' = L G™" and ^™" = L A™" one obtains 

(oji-n) a,7" = ^;"". (19) 

5™" is a new vector of Green's functions, each component 
T of which has only a single pole 



Amn.T 



(20) 



This allows us to apply the spectral theorem to each com- 
ponent separately. We introduce the vectors C™" = L C™" 
for the correlations, and 2?"™ = LD™" for the correction 
to the spectral theorem in case of a vanishing eigenvalue. 
Application of the spectral theorem 1 14 to the rth com- 
ponent of the single-pole Green's function of equation ( p^ 
then yields 



c- 



Artm.T 



1 



- + -(l-77)p'""■^ 



(21) 



where 



P" 






1 

2 

i lim i^^A^ 
2 tj^o u; ~ ll!. 



1 ^ r ^?TT.n,r 



(22) 



Here 6ro is the Kronecker symbol which ensures that 2?™"^"^ 
has a non-zero value only if r refers to the component with 
eigenvalue zero. 

Denoting L" as the left eigenvector corresponding to 
eigenvalue zero, we find 



V 



.17171,0 



^Ia. 



7nn.O 



1 



(Ap_ 



= - L'^ A™ 
2 C™") = 



, mn 

E 



i° c 



(23) 



Here we have used the relation between the commutator 
and anticommutator inhomogeneitics, A™"(k) = A™" + 
2C{^", and the regularity condition for the commutator 
Green's function for lo ^r Q 






(24) 



One sees explicitly, when inserting the left eigenvector of 
equation (HJ) belonging to eigenvalue zero, that 



E7-0 AOi.mn 1 



2El 



H-R'' ,H+H^-,2H^H' 



(25) 






are the regularity conditions of equation (17) of reference 
[||, see also Appendix B of this reference. 

The components of the correlation vector C™" for cUr ^ 
are of the form 



C" 



X 



7, = -l 



o/3wt — 1 



V 



.m7i,0 



(26) 



The original correlation vector C™" can be recovered from 
Qmn j-jy multiplying from the left with R; one obtains a 
compact expression by first defining a matrix C in terms 
of row-vectors C^ corresponding to the row-vectors of L : 



1 






L^ (r ^ 0) 



so that 



C™" = R C™" = R £ A™'_i + R P™" 



(27) 
(28) 

(29) 



The final equation determining the correlation vector is 
C™" = R £ A;'™_i + R" L" C™" . (30) 

The product R° L° is a projection operator onto the sub- 
space belonging to the eigenvector corresponding to to.,- = 

0, so that the term R*' L*' Qmn jg ^-j-^g projection of the 
correlation vector onto this subspace. This interpretation 
carries over to the iV— layer case, where, it will be seen, 
there is an iV— dimensional space corresponding to the 
zero eigenvalues. It is important to stress that this equa- 
tion is not complete but must be supported by the regular- 
ity conditions (Eq). Inserting the matrices R and L from 
equations (|lj) and (|l^) one sees that the z-component of 
this equation is exactly equation (27) of reference |g]. One 
could equivalently use the (-f ) or (— )-components of this 
equation, which can be proved to give the same results. 

In reference [0 we have investigated only spin S — 

1. In this case, it is sufficient to use the equations for 
[mn) — (01), (02) and (11). For general spin 5, all reg- 
ularity conditions with (to ^~ n) < 2S have to be taken 
into account. They form a set of linear equations which 
allow one to express all correlations ocurring in equation 
@) in terms of the moments ((S'^)^) {p ^ I, . . . ,2S + 1). 
This leads to 25* + 1 equations for the moments ((5*^)^), 
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Li Tij = -^^ Jij (H^ , Hi , 2H. 






{St) 

'h{st) 




^ J,, ( - H-{St) + HI{S-), H+iSt) - HIiSt), H~{St) ~ H+iST) 



(0,0,0) 



(41) 



which can be reduced to 25* equations by expressing the 
highest moment in terms of lower ones using the condi- 
tion Wj^j =-si^^ ~ Ms)=0. Note that only the first two 
equations have to be iterated because in the dispersion 
relation only (5^) and (S^S^) occur. For more details, see 
Appendix B. 



2.2 Multilayers 

Having established the formalism for the monolayer, it is 
now relatively easy to generalize to the multilayer case. 
For a ferromagnetic film with N layers the 3A^ equations 
of motion for the SN dimensional Green's function vector 
G'"" read 

(w 1 - r) G"" = A"" , (31) 

where 1 is the 3A^ x 3A^ unit matrix, and the Green's 
function and inhomogeneity vectors consist of N three- 
dimensional subvectors which are characterized by the 
layer indices i and j 










A Tiin 



'A. 



— jnn 



(32) 
The equations of motion arc then expressed in terms 
of these layer vectors, and 3x3 submatrices Fy of the 
2,N X 37V matrix T 



LUl- 



Tii ri2 

Tati rjV2 



TiAT 



G 



N] 



Aij(5ij 
A2j(5ij 



AnjSij 



(33) 

j = 1,..^N. After applying the decoupling procedures 
(0) and (^, the T matrix reduces to a band matrix with 
zeros in the Fy sub-matrices, when j > i -f 1 and j < 
i — \. The diagonal sub-matrices Tu are of size 3x3 and 
turn out to have the same structure as the matrix which 
characterizes the monolayer, see equation (pT[): 



(34) 



In particular one of the eigenvalues of Tu vanishes. The 
matrix elements of r^^ contain additional terms due to 
the exchange interaction between the atomic layers, the 
dipole coupling is contained in the field components Bf, 



Hf 





-H, 





-Hf 


H, 


hH- 


^H+ 

2^i 






see equation (13), 

H^ - Br + (sr) Ju {q ~ 7k) + jm+i( ^r+i) 

Hf = Hf + K2,, <?. ^ Z, + {Sf) Ju {q " 7k), 

Z, = Bf + A,+i (S^i) + J,,,_i (5f_i) + K2,r ^, , 

<^^ = 2{Sf)[l ^[S{S + 1) - {SfSf)]) , (35) 

and a = +, ~, z. The 3x3 non-diagonal sub-matrices Fy 
for j — HI are of the form 



~J,,{Sf) 

_ J,,(Sf| 







(36) 



We now demonstrate that there is a left eigenvector of F 
corresponding to eigenvalue zero with the structure 



LO = (0,...,0,LO,0,...,0), 
where 

tO _ /rO 7-0 7-0 



(37) 



1 



2El 



H^Hf, H+Hf, 2HfHl 



(38) 

This is immediately clear for the diagonal elements, be- 
cause they have the same structure as the monolayer ma- 
trix, equation (O). To prove this also for the non-diagonal 
matrix elements one needs the regularity condition ( |2q ) for 
layer i: 

Lt Af '"" = , (39) 



E 



for m = 0, n 



A' 



1. With Af'°^ = 2(5f}, Af'"^ = 0, and 



-{S^ ) we obtain 
{Sf) = 



Hf 



{SD 



(40) 



With this regularity condition, we complete the proof that 
L*' is a left eigenvector of F with eigenvalue zero since we 
obtain for the non-diagonal elements 

see equation (pl| ) above. 

Hence, N out of the 3A^ eigenvalues of the multilayer ma- 
trix F must be zero. 

Apart from dimension, the equations determining the 
correlation functions have the same form as for the mono- 
layer case: 



R£A;™_i 



■p T r~imn 



(42) 
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The matrices R and L have to be constructed from the 
right and left eigenvectors corresponding to non-zero eigen- 
values as before, whereas the matrices R" and L° are con- 
structed from the N eigenvectors with eigenvalues zero. 

In order to compute the matrix T when iterating equa- 
tions (^) , one has to solve the linear system of equations 
(1), 



Z,{Sf) ~Ju+i {St,,) (5f ) 



Jv.-i {St^) {S!) 



= B^ (S^) 



z = l,...,iV, 



(43) 



in each iteration step. For general spin S, one can express 
all higher correlations occuring in equation (B3) in terms of 
the moments of {{S^)'^^'^^) for layer i by using all regular- 
ity conditions with m + n < 2S. Again the largest moment 
can be expressed by lower ones using JJj^j [Sf — Ms) = 0. 



2.3 The effective anisotropy 

If, in addition to the orientation angle, one is interested 
in the effective (temperature-dependent) anisotropy coeffi- 
cient K2{T), a quantity which is accessible in experiment, 
one needs a working expression for the free energy. For 
a derivation of this expression, which is also used by ex- 
perimentalists to extract K2{T), we refer to the books of 
Landau and Lifschitz ||l5[ and of Vonsovskii ijl^l . To lowest 
order the free energy reads 

N 

F{T) = J2 F^{T) 
1=1 

HT) = ~\Y1 -^a S, • Si - K2,m cos2 e,; - B • S, 
<i> 

+ ^Ef?(s^-S<-3(S,.r,0(S,-r,0).(44) 



I 



As in reference B , the temperature-dependent anisotr- 
opy K2.i{T) for each layer i is calculated non-perturbatively 
by minimizing the free energy with respect to the layer- 
dependent reorientation angle 0^. From the condition 
dF{T)/de,{eQ,^) = 0, we find with B = (S=^, 0, B^) 



Nh{T) 



is the exchange interaction between neighboring layers, 
and 

^-^=E5^.^(!';M^'^-^'' (46) 



2sm6'o,i cosfeio,; 
cos 6>o,j (^S"" + Jj,j+i Mj+i(T) sineci+i 

+J,,,_i Af,_i(T) sineo,»-i +Tr) 

- sinOa^iyB'- + J,,,+i M,+i{T) cos6>o,i-Hi 

+J,,,_iA/^_i(T) cos 00,^-1 -2T^ 



(45) 



where the magnetization M,{T) = ^(S'f )2 + (5f )2, and 
the equihbrium polar angle ©o.i = arctan((5f )/(S'f )) are 
determined from the magnetization components {Sf) and 
{Sf) calculated from the Green's function method. Ji.i±i 



cosOq,. 



TI' ■'I being the dipole lattice sum occuring in equation 
(Esh. For a single layer A^ = 1, we obtain equation (31) of 
reference ||, with an additional term (3/2) gn Mf (T) T° 
due to the dipole coupling. The total effective anisotropy 
K2 (T) of the thin film is given by 



N 



X2(T)=^if2,.(T) 



(47) 



This procedure for the determination of the effective anis- 
otropies K2,i{T) is non-perturbative in the sense that the 
magnetization and the orientation angle in equation ( |4q ) 
are calculated from the full Hamiltonian, in contrast to a 
thermodynamic perturbation theory where one splits the 
Hamiltonian into two terms, e.g. |l[. In the numerical cal- 
culations, we will normalise the anisotropy coefficient in 
the Hamiltonian to K2,i/S{S — 1/2), in order to guarantee 
that K24T)/K24T = 0) = 1 for T ^ 0. 

We are aware of the fact that using only K2 in the 
Hamiltonian can lead to an effective KiiT) Q], which, al- 
though it turns out to be very small in an analysis within 
mean field theory, has some effect on the nature of the 
phase transition and on the phase diagram. We do not try 
to extract a corresponding K4 term here, because we do 
not calculate the order of the reorientation phase transi- 
tion or a phase diagram in the present paper. 



3 Results 

In this section we show results of the calculations de- 
scribed above. First we discuss the effect of the dipole 
coupling on the reorientation of the magnetization for a 
monolayer with spin S = 1. Secondly, we discuss the case 
of a single layer with spins S > 1. Thirdly, we treat ferro- 
magnetic films consisting of A^ layers. 



3.1 The effect of the dipole coupling on the magnetic 
reorientation of a monolayer 

Since in reference [gj the dipole coupling has not been 
taken into account explicitly, we investigate in this subsec- 
tion the action of this interaction on the magnetic reorien- 
tation in the case of a monolayer with spin quantum num- 
ber S = 1. As discussed above, the exchange coupling is 
treated by RPA, the single-ion anisotropy terms according 
to the Anderson-Callen decoupling, and the dipole cou- 
pling is considered within the simplified (non-dispersive 
mean field) approximation described in Appendix A. We 
use the parameters J = 100, K2 = 1 chosen in reference 
[||. The dipole couphng strength is set equal to gn = g = 
0.018 or g = 0.066, which refers to the cases of Ni or Co 
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Fig. 1. The components of the magnetization (S^) (solid lines) 
and {S^) (dashed lines) of a spin S = 1 monolayer are shown 
as functions of the external magnetic field in x-direction, B^ , 
without [g = 0) and with dipole coupling [g = 0.018, the value 
estimated for Ni), and J = 100 and K2 = 1. The dipole cou- 
pling renormalizes the external magnetic field. Two different 
reduced temperatures (T jTc = 0.5 and 0.8) are considered, 
where Tc is the Curie temperature for perpendicular mag- 
netization. The (S^) components are plotted only up to the 
temperature where (5*^) -^ 0. 
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Fig. 2. Effective anisotropy K2{T)/K2{0) and equilibrium 
reorientation angle 0o of a spin 5 = 1 monolayer are shown 
as functions of the temperature T without (g = 0) and with 
dipole coupling {g = 0.018 for Ni and g — 0.066 for Co). A 
magnetic field B^ = 0.1 is applied, and J = 100 and K2 = 1- 



by calculating the relative strength of J/g, where J is es- 
timated from the corresponding bulk Curie temperatures. 
The external magnetic field is directed along the x-axis, 
B = (5^,0,0). 

In Figure 1 we plot the components of the magnetiza- 
tion (S*^) and (5'=") without {g = 0) and with (.9 = 0.018) 
dipole coupling for the temperatures T/Tq — 0.5 and 
T /Tc — 0.8 as a function of the transverse field B^ . Here 



Fig. 3. Normalized magnetizations {S^)/S and {S'^)/S for 
a monolayer as functions of the temperature for all integral 
and half-integral values of the spin between S = 1 and 5 = 6 
calculated with the Green's function theory. The reorientation 
temperature Tf depends slightly on S. The inset shows corre- 
sponding results where all interactions are treated with MFT 
for spins S = 1,2,7/2, and 11/2. In this case, the reorienta- 
tion temperature T^^'^ does not depend on 5*. All parameters, 
B^ = 0.1, J = 100, K2 = 1, and g = 0.018, are scaled as dis- 
cussed in the text. 



Tc is the Curie temperature calculated for a perpendic- 
ular magnetization. As expected from equations (|l2|,|l3|), 
the dipole coupling diminishes the action of the uniax- 
ial anisotropy and the magnetic field in the z-direction, 
leading to a reduction of (S^), and enhances the action of 
the transverse components; consequently, {S'^) increases. 
Therefore, the reorientation field B^, at which (S*^) van- 
ishes, becomes smaller for increasing dipole coupling 
strength. 

In Figure 2 we plot the equilibrium reorientation an- 
gle 6)o(T) = arctan((5^)/(5'^)) for S^ = 0.1 as a func- 
tion of the temperature T. The following dipole coupling 
strengths are considered: g = 0, g — 0.018 (estimated 
for Ni), and g = 0.066 (estimated for Co). With increas- 
ing dipole coupling strength the reorientation temperature 
Tfj, at which {S") vanishes (6*0 = 90°), decreases. The cor- 
responding effective (temperature-dependent) anisotropy 
K2(T)/K2{0) as obtained from equation (|45|), is also shown. 
Since the results shown in Figure 2 are obtained for a finite 
magnetic field B^ = 0.1, (S^) remains finite for T > Tr, 
and K2{T) does not vanish completely at Tr. For the cou- 
pling constants under consideration the overall behaviour 
of K2{T)/K2{0) does not change for different g. 



3.2 The monolayer for S > 1 

In this subsection we investigate the effect of different spin 
quantum numbers S on the magnetic reorientation. We 
consider a monolayer with the interaction parameters used 
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Fig. 4. Results of a Green's function calculation for the 
normahzed moments ((5'^)")/5'" for a monolayer and spin 
S = 11/2 for n — 1,2, ...,11. Tr is the reorientation tem- 
perature. 



for the results of Figure 1. In order to compare the results 
for different S we have scaled these parameters in the fol- 
lowing way: J ^ J/S{S+1), B -> B/S, g -^ g/S{S + l), 
and K2 = K2/ S{S — 1/2). The scaling of K2 guarantees 
the property \imT-,oK2{T)/ K2{Q) = 1 for B^ -^ 0. 

In Figure 3 we display the normalized magnetizations 
{S^)/S as functions of the temperature T for all inte- 
gral and half-integral spins ranging between 5=1 and 
S = 6. The reorientation temperature, Tf{{S), becomes 
smaller with increasing S but at a rapidly decreasing rate 
as 5* increases. The small external magnetic field in x- 
direction (B^ = 0.1) induces a finite x-component of the 
magnetization {S^)/S for T = 0, which increases slightly 
with increasing temperature. These results are compared 
with results of calculations where a mean field approx- 
imation (MFT) is performed for all interactions [p[, see 
the inset of Figure 3. Within this aproximation a more 
pronounced spin dependence of the magnetization curves 
is observed. On the other hand, due to the scaling of the 
coupling parameters the mean field reorientation tempera- 
tures T^^'^{S) are independent of the spin value S. Note, 
however, that T^^^ is more than a factor of two larger 
than the reorientation temperature as calculated from the 
Green's function theory - this is due to missing correla- 
tions in MFT. 



In Appendix B it is shown that the correlations ocur- 
ring in the equations of motion can be expressed by higher 
moments of the magnetization, which are determined by 
the regularity conditions. As an example, we present in 
Figure 4 results of the spin wave calculation for the nor- 
malized moments {{S^)^)/S^ for spin S = 11/2 and for 
n = 1,2,.. .,11 as functions of the temperature. The odd 
moments approach zero for T —>■ Tn, whereas the even 
moments approach a finite value for T — > Tg. as expected 
physically. For n = 2 one obtains for example ((S'^)^)/5^ — > 
S{S + l)/iS'^. 



In Figure 5 the equilibrium reorientation angles ©o, 
and in Figure 6 the corresponding effective anisotropy co- 
efficients K2{T)/K2{0) are shown as functions of the tem- 
perature for spin quantum numbers ranging from S — 1 
to S = 6, as calculated from the Green's function method. 
The temperature dependence of both these quantites does 
not vary markedly with S. In the insets of the figures are 
the corresponding results for 0q and K2{T) / K2{0) as ob- 
tained from a MFT for all interactions. As already seen 
for the magnetizations, a more pronounced spin depen- 
dence of K2{T)/K2{0) is observed here also for MFT. Ev- 
idently, the spin quantum number 5* has a larger influence 
on single-spin excitations (MFT) than on collective mag- 
netic excitations (spin waves). 

If scaled coupling constants are used, only a weak de- 
pendence of the spin quantum number S on the magnetic 
quantities such as the magnetization, the reorientation an- 
gle, and the effective anisotropy is observed within the 
Green's function method. Thus, one may perform calcu- 
lations with a low spin, for which a considerably smaller 
system of equations has to be solved self-consistently. Re- 
sults for higher spins can then be obtained by scaling. This 
is less justified within the MFT approximation. 

We stress a result already obtained in reference 12] for 
spin S — 1, namely that the effective anisotropics K2{T) 
as calculated within the RPA and within MFT have dif- 
ferent temperature behaviours, particularly at low tem- 
peratures. In this temperature regime, the MFT exhibits 
an exponential behaviour, and RPA an almost linear be- 
haviour of K2{T). Consequently, the use of MFT would 
lead to a considerably smaller value for -^'2(0) than that 
obtained with RPA, when the observed values of K2{T) 
(measured typically at T/Tc — 0.7, e.g. Q) are extrap- 
olated to T = 0. Note that only at T = are anisotropy 
coefficients available from ab-initio calculations, e.g. refer- 
ence Q. These theoretical values can only be compared 
with experiment by extrapolating measurements at finite 
temperatures down to zero with the help of a theoreti- 
cal model. Because of the drawbacks of MFT we propose 
performing this extrapolation with the results from the 
Green's function method. 



3.3 Ferromagnetic films with N layers 

In this subsection we demonstrate that the formalism de- 
veloped in Section 2 can be applied to the case of many 
layers. We study the magnetization, the reorientation an- 
gle, and the effective anisotropy as functions of the film 
thickness (characterized by the number of atomic layers, 
N). In all the examples shown in this subsection, the 
dipole coupling is included within the simplified mean field 
treatment as discussed above. 

As examples we treat simple cubic films consisting of 
N layers with spin S — 1, using the same coupling parame- 
ters for all film layers (homogeneous film): J = Jih = 100, 
K2,k = K2 = 1, gik = g = 0.018. These couplings are 
scaled as in Section 3.2 : J ^ J/S{S + l), g -^ g/S{S+l), 
K2 -^ K2/S{S - 1/2), and B -^ B/S. All quantities 



p. Frobrich et al.: Many-body Green's function theory for the magnetic reorientation 




Fig. 5. Equilibrium reorientation angle 6*0 as a function of 
the temperature for the systems of Figure 3 calculated with 
the Green's function theory. The inset shows the corresponding 
results when applying MFT to all interactions. 
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Fig. 6. Effective anisotropy K2{T)/K2{0) calculated with the 
Green's function method as a function of the temperature for 
the systems of Figure 3. The inset shows the corresponding 
MFT results. 



are calculated for a small external magnetic field B = 
(0.1,0,0). 

In Figure 7 we show the sublayer magnetizations (5*^^), 
i = l,...,iV, as functions of the temperature for film 
thicknesses ranging between A^ = 1 and A^ = 19 lay- 
ers. As expected, and also seen in MFT calculations [^, 
m , for a homogeneous film, the magnetization of the sur- 
face layers is smaller than those of the interior layers be- 
cause of the smaller coordination number for the surface. 
Also, the magnetizations of the ith and the {N — i + l)th 
layer are the same (twofold symmetry). For the parame- 
ters under consideration the reorientation temperature Tn 
is close to the Curie temperature for perpendicular mag- 
netization. Therefore, we see in Figure 7 that the value 
for the reorientation temperature Tn (defined as the tem- 



Fig. 7. Sublayer magnetizations (Sf) as functions of the tem- 
perature for thin ferromagnetic films with A'' layers and spin 
S = 1. The reorientation temperatures T^ for the different 
films can be read off from the curve in the N — T plane, 
where (Sf) — 0. The same parameters are used for all lay- 
ers: B^ = 0.1, J = 100, K2 = 1, and g = 0.018. They are 
scaled as described in the text. 



perature where (5f ) = 0) exhibits a saturation behaviour 
as a function of the film thickness. Whereas there is a 
steep rise from Tii{N = 1) = 58.55 for the monolayer 
to Tb.{N = 3) = 103.08 for the trilayer, there is only 
a small difference between Tfi{N = 17) = 135.11 and 
Tr{N = 19) = 135.67. 

Because in experiment only the average orientation of 
the thin film magnetization is measured one has to cal- 
culate this quantity from the model. In Figure 8 we show 
the average equilibrium reorientation angles 6*0 (iV, T) of 
thin films with different thicknesses N as functions of the 
temperature, where 



OoiN,T) = arctan 






(48) 



An alternative to calculating the average orientation an- 
gle is first to calculate the angles for each layer, 0o,i = 
arctan((S'f)/(S'f)), and then to average over the angles. 
The difference between both procedures turns out to be 
tiny. 

In Figure 9 the average effective anisotropics K2 {N, T) 
/K2{N,0) of thin films with different thicknesses N, cal- 
culated from equation (47), are shown as functions of the 
temperature. With increasing film thickness the action of 
the effective anisotropy extends to higher temperatures. 
In the inset we also show for iV = 1 and TV = 19 that 
the dependence of the anisotropics on the temperature 
is somewhat different for different layer thicknesses if one 
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Fig. 8. The average equilibrium reorientation angle ©o as 
a function of the temperature and film thickness. A^ is the 
number of layers in each film and T^ are the reorientation 
temperatures (6>o = 90°). 
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Fig. 9. Average effective anisotropy K-^.i^') j K^iS*) as a func- 
tion of the temperature and film thickness A'^. The inset demon- 
strates the different functional dependence of K2{T) / K2{0) on 
the temperature for layers with A^ = 1 and A*' = 19 if the 
temperature T is scaled with respect to the reorientation tem- 
perature Th- 



scales the temperature to the respective reorientation tem- 
peratures: T/Tii{N). The curvature of the curve for the 
thick film (TV = 19) is more pronounced as that for the 
monolayer. 



4 Discussion and conclusions 

In thepresent paper we have extended the model of refer- 
ence [g] in various respects and we have developed a suc- 
cinct formulation of the final equations. This was made 
possible by utilizing the eigenvectors as well as the eigen- 
values of the matrix governing the equations of motion 
for the set of Green's functions, which has to be intro- 
duced when the calculation of several non- vanishing com- 



ponents of the magnetization is required. This new pro- 
cedure, which for the monolayer and spin S" = 1 is fully 
equivalent to our earlier treatment |g], provides a prac- 
ticable way of extending the Green's function spin wave 
theory to the reorientation of the magnetization of ferro- 
magnetic films consisting of many layers and for general 
spin S. 

We have applied the new method to the monolayer case 
with spins S > 1. We have found that the spin dependence 
of the magnetizations and the anisotropics as functions of 
the temperature is considerably less pronounced in RPA 
than in MFT if a proper scaling of the parameters of the 
model is applied . The corresponding curves saturate much 
more quickly in RPA than in MFT with increasing spin 
quantum number S. The temperature for complete reori- 
entation Tr, on the other hand, does not change with S 
in MFT, whereas there remains a slight spin dependence 
in RPA. 

For the monolayer with spin 5 = 1, we have investi- 
gated in detail the influence of the dipole coupling on the 
reorientation problem. Because, for realistic dipole cou- 
pling strengths, we found no big differences in treating 
the dipole coupling with its RPA or (non-dispersive) MFT 
approximations , cf. Appendix A, we chose to include the 
dipole coupling by means of the latter, which is relatively 
simple to handle and requires only a renormalization of 
the external magnetic field. 

We emphasize that only by using our new method have 
we been able to treat the magnetic reorientation within a 
Green's function approach for films with several layers. We 
have studied the field-induced magnetic reorientation and 
the effective anisotropy as functions of the film thickness 
and temperature for spin 5=1 films. Investigations of 
films with S > 1 present no problem; they arc only more 
time consuming. 

In the present paper, we have only studied homoge- 
neous films, but the method can also treat inhomoge- 
neous films or multilayers by using layer-dependent cou- 
pling constants and magnetic moments. The magnetic re- 
orientation could then be calculated for thin film or mul- 
tilayer systems investigated experimentally. This will be 
pursued in forthcoming studies. 

A few words concerning the applications of the present 
model are in order. A prerequisite is that the investigated 
systems can be modelled in a reasonable way by a local 
spin model of Heisenberg type. Moreover, it is required 
that higher order anisotropics are not important, a condi- 
tion which is often not fulfilled. 

Most results of the present paper are obtained for a 
small transverse field which initializes the reorientation. 
It is certainly of interest to calculate also phase diagrams 
in which the magnetic field is varied, and to compare with 
results obtained in reference |2(]|l within a schematic model 



and in reference |21 on the basis of mean field theory. 

The possibility of treating spins 5 > 1 would also al- 
low the treatment of the fourth-order uniaxial anisotropy 
K4 as well as the quartic in-plane anisotropy. This, how- 
ever, requires a proper decoupling procedure for the cor- 
responding terms in the Green's function theory which we 



p. Frobrich et al.: Many-body Green's function theory for the magnetic reorientation 



11 



do not have available at the moment. A phase diagram 
in the K2 — K^ plane would then show the region of sta- 
ble magnetization directions and would help in exploring 
the locations of temperature or film thickness driven mag- 
netic reorientations. Also, one could determine whether 
the magnetic reorientation happens continuously or dis- 
continuously. 



Appendix A: The treatment of the dipole cou- 
pling 

In this Appendix we show how the long-range magnetic 
dipole coupling can be considered within the RPA treat- 
ment of the magnetic reorientation. Furthermore, we show 
that a simplified treatment for interaction strengths small 
compared to the exchange interaction leads to a satisfac- 
tory description of the magnetic properties. For simplicity 
we consider only the case of a single layer {N = 1). 

As it should be, the magnetic dipole interaction leads 
to an additional dispersion in the magnon dispersion re- 
lation a-'(k). Applying the generalized Tyablikov (RPA) 
decoupling, equation (0), to the dipole interaction in the 
Green's function equations of motion (||), ( ( [S'f , W^'p°'°] - ; 
{Sj)"^{S~)^)), one obtains the following additional terms 
on the left side of equations ( |lO| ) 



-T^ -Tu -T^ 




(49) 



Here a 2D Fourier transformation has been applied, and 

^k = ^ 5 {S'){T^, T^2 + 2* n,) (50) 

Ti = g (5+) [n, + T^2 + \n, + iT°2 



where 



-i{S+){T^^^T^2-2tT^^, 






hn 






V: 



2 N5/2 



exp{ikxXi) ex.p{ikyym) (51) 



are oscillating lattice sums, which can be evaluated with 
Ewald summation techniques as outlined for instance in 
reference Q. 

As mentioned in the main body of the paper, the RPA 
treatment of the magnetic dipole coupling complicates the 
calculation of the magnetization considerably because of 
the presence of complex terms and dispersive (k— depen- 
dent) terms. Thus, as an approximation we neglect now 
the dispersive parts in the equations of motion (^) coming 
from the dipole coupling, and retain the non-dispersive 



terms only. This corresponds to a mean field treatment of 
the dipole coupling. Then equations (pO|) reduce to 



T.^ 



9{S')iK + T^2) 



T-=0 



(52) 



iS+KK+To': 



02) 



Ti^-~i{S-){n+T^2)- 

This simplification allows the dipole coupling to be taken 
into account by a renormalization of the external magnetic 
field, and leads to equations ( p^ ) of Section 2. 

By comparing to RPA, we shall show now that this ap- 
proximation leads to satisfactory results for small dipole 
coupling strengths as found for instance in ferromagnetic 
3d- transition metal thin films. Because the general RPA 
treatment of the dipole coupling turns out to be fairly 
complicated, we consider only two limiting cases which 
are manageable, a perpendicular and an in-plane magne- 
tization. An external magnetic field is not considered but 
can be easily added. 

For spin S' = 1, one needs the Green's functions G^ ■'"*" 

= ((5'=^;(5'|)'"(S'-)")) for n = 1 and for m = and 
m — 1. In case of the perpendicular magnetization, use 
of the RPA decoupling for the exchange interaction and 
the dipole coupling, and the Anderson-Callen decoupling 
for the single-ion anisotropy (<? — {S^){S^S^) for 5 = 1) 
leads to the following equations of motion (the 2-axis is 
directed perpendicular to the plane) 



-b* cu + aj \G-'"'\k,J) 






(53) 



with 



« = {SI (j (q 7k) - 9{t^o + ^02 + It^o + It^2) 



-K2 {S^S'' 



b = 



l9{S')[T^,-T^2 + '^^T^i 



(54) 



Solving these equations for the Green's functions and ap- 
plying the spectral theorem we obtain the following cor- 
relation functions 



+ L{aAtr^-bA-_r'^cot\,^^' 

{{STS-S-) = ]-Azr^ 



i'i)^ 



1 / /I -.ml 

2^r^-i 



b* Atr^ ) coth 



(f), (») 



where the dispersion relation is 



(56) 
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Using now for S" = 1 the identities S-S+ ^2-S^- {S^'f 
and {S^)'^ — S^, the inhomogcneities are given by 



A 



+m 



1+41 






2(5'") , AZ\ 
i-.ii 



(s-^ 



(57) 



For m — we obtain 

((5^)2) =2- (5") (1 + 201), 

and for m = 1 



{S-S-) = 2{S^)<I>;, 
(58) 



(3((5")^)-(5")-2)0i-(^-^- 



(P2 , 

(59) 



where sums (integrals) over the first Brillouin zone have 
to be performed for the quantities 



N '^\2€ ^ 2 ' 



62 — — > — coth — . 
^ N ^2t ^ 2 ' 

k 

From these equations we find the expectation values 

1 + 201 



{sn = 



l + 30i+302+|02|2 



((5")2) = 2-(5")(l + 20i) 



(60) 

(61) 
(62) 



These equations are solved self-consistently for (S^) and 
((5^)2). 

In Figure 10 wc show these expectation values with full 
RPA for the dipole coupling {g = 0.018 estimated for Ni) 
for a monolayer using the parameters J — 100, K2 = 1, 
and for g — (no dipole coupling). For the perpendicular 
magnetization, the {S^ S^) correlation is very small, and, 
within the line thickness of the curves, there is no differ- 
ence between the full RPA and a RPA where one puts 
(S^S^) = (6 = 0). Also shown are the corresponding 
results for the dipole coupling considered within the (non- 
dispersive) mean field treatment. In this case, not only 
6 = but also the k-dependent terms in equation ( p4[ ) 
connected with the dipole coupling are neglected, i.e. 

« = (S^) {j {q 7k) - 9 (T°„ + T,%) K^ {S^ S^)) , 
5 = 0. (63) 

As can be seen from Figure 10, the difference between 
a RPA and a MFT treatment of the dipole coupling is 
small. 

The situation is somewhat different for an in-plane 
magnetization. In this case the single-ion anisotropy is not 
active (^ = 0), and the dipole term is, in accordance with 
the Mermin- Wagner theorem, the only term which induces 
a finite magnetization. In an RPA treatment of the dipole 
coupling the magnetization in the plane is not isotropic 
and one has to introduce an in-plane angle with respect 
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Fig. 10. Expectation value (5*^) for a monolayer with perpen- 
dicular magnetization and for spin 5 = 1 as a function of the 
temperature. The dipole c oup ling [g = 0.018 for Ni) is treated 
with full RPA, equation (pj), and its mean field approxima- 
tion, equation (p3|). For comparison we show also the result 
with vanishing dipole coupling (g = 0). We have used J = 100, 
K2 = 1, and a vanishing external field B = 0. 
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Fig. 11. In-plane magnetization (5^) for a monolayer with 
spin S = 1 as function of the temperature. The dipole cou- 
pling {g — 0.018) is treated with full RPA, equation (|6J), with 
RPA putting (S^S^) = 0, i.e. 6 = in equation (pq), and 
with MFT, equation (pq). The {S~ S~)rpa correlation is also 
shown. 



to a main axis of the square lattice. By putting this angle 
equal to zero one obtains for the quantities a and b, cf. 
equation (|5S 



(^^)(^('Z-7k) 



{T^o 



n 



n 



2n 



b=^9 (5") T, 



(64) 
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Using the mean field approximation for the dipole cou- 
pHng by neglecting the corresponding k-dependent terms 
one obtains for the in-plane magnetization 



(5^>(^('?-7k) + f (tO„ + T°2)) 







(65) 



Within the formalism above, the quantization z-axis is 
now directed in-plane, which corresponds to the s-direction 
of our original reorientation problem. The corresponding 
results are shown in Figure 1 1 . Here the difference between 
the full RPA and MFT is somewhat larger, because the 
expectation value {S^ S^) is not as small as in the case for 
the perpendicular magnetization. In particular the Curie 
temperature is appreciably lower than in MFT. This, how- 
ever, is the most extreme case, and is an upper limit for 
the error one makes if the dipole coupling is treated by 
MFT for the reorientation problem. As also shown in Fig- 
ure 11 the MFT result is close to the RPA result only 
when the {S~ S~) correlation is neglected. 

In summary, the error is small for the perpendicular 
orientation and will increase with increasing polar angle, 
but is expected to stay below the error for the in-plane 
case when replacing the full RPA treatment of the dipole 
coupling by the simplified, non-dispersive MFT. With this 
discussion in mind we use for the reorientation problem of 
the present paper the mean field approximation for the 
dipole coupling. We do so also because its full RPA treat- 
ment is quite complicated and time-consuming. If the ratio 
between the dipole coupling strength and the exchange 
coupling becomes larger, as expected e.g. for rare earth 
ferromagnets, the error will increase and the simplified 
treatment is less justified. In this case one is faced with 
the more complicated full RPA treatment of the dipole 
coupling. 

Finally, we mention that the magnetic dipole coupling 
can be taken into account also via its mean field approxi- 
mation when calculating the magnetic reorientation with 
the Schwinger-Boson theory [|lO| . However, the proper 
inclusion of the long-range character of this interaction, 
which leads to additional momentum-dependent terms in 
the magnon dispersion relation, is practically impossible 
within this method. 



Appendix B: Treating S > 1 

In this Appendix we show how the regularity conditions, 
which have been calculated in [|| for 5 = 1, can be de- 
duced for general spin quantum numbers S. The regularity 
condition (G6h for m — 0,n — l yields 



H-^ Z ' 



(66) 



and can therefore be written for general ?7i, n in the form 



iZA'' 



A 



+ ,77in rj — 



a:j:'""b+ 



(67) 



The z-component of equation (pG) , from which the corre- 
lations have to be calculated, reads then explicitly 

~ ^^^ (|(5")"(5-)"5^) - ^ ^(5")™(5-)"5+) 
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(68) 
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We express all correlation functions occuring in this equa- 
tion in a standard form where all powers of S^ are written 
to the left of the powers of S^ : 



C{m,n)::^ ((5^)"(5'-)") 



(69) 



Then, with the relations [5^(S•-)"]_ = -n(5-)" and 
S-S+ = S{S +1)-S^- (5*^)2, we find that 

{{S'nS-yS') =nC(m,n) + C(m + l,n) , 
{iSTiS-y'S+) = {S{S +l)-n{n- 1)) C{m, n - 1) 

~{2n~l)C{m+l,n- 1) - C(m + 2, n - 1) , 
{{STiSTS') = Cim,n+l). (70) 

The commutators can also be expressed in terms of the 
C{jn, n) using the binomial series 
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n 
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(71) 



Now by putting equation (|7l]) into equation (|6^) the reg- 
ularity conditions for all m and n can be written in terms 
of correlations defined in the standard form: 



2 Z n C{m, n) 



B 



{-\yC{m-i,n-\) 



m+l 



5(5 + l)E 

i=l 

+ Y.{"t^){-^r'C{m + 2-^,n-l) 

+ {2n + m) C{m + 1, n - 1) -|- n (n - 1) C{m, n - 1) 

+B^il{'l:) C(m-z,n-fl). (72) 
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For a given spin S*, this set of Unear equations for the 
correlations has to be solved for all m + rt < 25+1. 
The solutions have to be put via equations (^) together 
with ( [n]) into equations (|6q), thus leading to a set of 25 
equations for the moments {{S^Y) (p = 1, . . . , 25), which 
have to be solved self-consistently. The highest moment 
{{S^Y^^^) has been eliminated in favour of the lower ones 
through the relation Hji/s (5"' - Ms) = 0. 
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